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Abstract 

We propose a matrix-model derivation of the scaling exponents of the critical and 
tricritical q-states Potts model coupled to gravity on a sphere. In close analogy with 
the 0(n) model, we reduce the determination of the one-loop-to- vacuum expectation 
to the resolution of algebraic equations; and find the explicit scaling law for the case 
q=3. 
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1 Matrix realization of the Potts model 



Matrix models have proved to be an efficient means to deal with the statistical properties 
of matter coupled to two-dimensional gravity: instead of performing functional averages 
on the metric and embedding fields, we consider a matrix-valued field theory whose dia- 
grammatics corresponds to the enumeration of simplicial complexes (discretized surfaces) 
dressed with other (i.e. non-geometrical) degrees of freedom (e.g. Ising, Potts spins, a 
scalar field...). 

Here, we shall examine 0: 



with 0i. ..g q N x N hermitean matrices, and (ij) denoting any pair of different indices. 
Expanding in powers of g and performing the gaussian integral, we obtain various 3 
fat graphs (we represent any propagator with a double-line and corresponding row/column 
indices attached, and any vertex with a fat Y with the index structure of the trace of a 
third power <p a b ( t ) bc ( t ) ca) where each vertex has been attached a number i = 1, . . . q ( which 
matrix has been chosen to produce that vertex?). We can as well consider dual graphs: 
various gluings of colored triangles in an orientable surface without boundary. As for the 
statistical weight, we get a contribution from the vertices, recognized as a cosmological 
constant (energy~area = number of triangles); and from the propagators,e.g.: 



where the energy of a satisfied link (neighbouring triangles with the same color) re- 
normalizes the cosmological constant (there is 1.5 link per unit area) and the relative 
energy of frustrated bonds produces the required statistical weight of the Potts model; 
the temperature is truely positive, provided that m 2 > q. For any m, the average size 
of the surfaces will grow with g, and there is a critical value g c {fn) near which the ther- 
modynamic regime is achieved; if we also tune m, the temperature, the spin-correlation 
length will also diverge, and we shall obtain the critical regime. 

So far, diagrams of all genera appear; but they are weighted by N x=Euler charactensUc 
as is easily seen: a factor of N accompanies any vertex; 1/N for each propagator; and 
a choice of a free index (a = 1,...N) for every loop in the 3 fat graph. Taking the 
logarithm, we obtain connected diagrams only, and sphere-like diagrams are extracted in 
the N — > oo limit. 

In all solved matrix models, a tremendous simplification appears in that limit: vari- 
ables are changed from to Q, A according to = fi + diag 1 ^(A)f2 and when the integral 
over the diagonalizing unitary matrices is preformed, one is left with an integral over the 
N eigenvalues per original matrix with an effective energy O(NTr) = 0(N 2 ), that is: a 
semi-classical problem (note that the ground state corresponds to eigenvalues of order 1). 
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Indeed, the integral over Q is trivial for a one-matrix model where the energy is 
invariant under conjugation. For terms coupling two matrices 

Tr0i0 2 = ™+diag(Ai)fii^diag(A 2 )fi 2 

the required integral involves the relative angle Q1Q2 an d is exactly known |J. However, 
when the model couples 3 or more matrices in a cyclic way {4>i4>2, <p2<j>3, </>30i) the relative 
angles are no more independent, the semi-classical reduction is no longer possible and we 
do not know how to solve such models - that would correspond to matter with c > 1 (e.g. 
a 2-dimensionnal lattice of matrices to produce surfaces embedded in 2 dimensions). 

This problem is only apparent in the case under consideration: the Potts model on a 
flat regular lattice is known to have central charge less than one, and the cyclic structure 
of ([1]) can be disentangled by introducing an auxiliary gaussian matrix symmetrically 
coupled to the Potts matrices [|l[ [j 

Z = J dXdtexp-NTr^Vtyi) +X 2 /2 -^fax] (3) 

i i 

where V{<j>) = ^(f) 2 + § 3 . 

Now the relative angles Slx/fc are independent and can be integrated out: to solve 
the problem for planar random graphs, we just have to determine the distribution of the 
eigenvalues of the auxiliary and Potts matrices in the lowest-energy configuration. We 
shall note p the density of eigenvalues for X, that is p(x) dx = the fraction of eigenvalues 
in [x, x + dx] ; a will correspond to the Potts matrices (we assume no symmetry breaking, 
so that all q matrices have the same spectrum at equilibrium). The associated resolvents 
are / and g: 

(4) 

and correspondingly for g; they are holomorphic functions on the complex plane cut along 
the supports of their respective densities (finite intervals). The densities do depend on m 
and g, and scaling laws exist between the critical exponents attached to them and critical 
indices of the spin system: for instance, the string susceptibility 7,^ which governs the 
area-dependence of the two-point function is known through the exponent with which the 
density of eigenvalues vanishes near the edge of its support exactly at the critical point 
a(y) ~ y 5 

Istr = "(5 " 1) (5) 

We shall find the densities in the critical regime. 

When we change from hermitean to polar variables, a Jacobian A(A) 2 appears: the 
square of the Vandermonde determinant of the eigenvalues. This leads to (repulsive) 
interactions between the eigenvalues, to which we shall add the effect of the potential 
dressed by the interactions between matrices: if we set 

/(X) = f d<Pe- NTl ^- x ^ = I( Xl , ...,x N ) (6) 



the dressed potential seen by the x's is NJ2 x t — <l log (H x )) an d the classical equilibrium 
equation reads 

2 1 

— J2 — Xi + qw(xi) = (7) 

IV Xi Xj 

where w{xj) = log(J) = (fin), the average being taken for a matrix (f> that fluctuates 
in its own cubic potential under the influence of a diagonal matrix X with spectrum p. 
So, we have 

2Re f(x) — x + qw(x) = , x G supp p (8) 

and if we introduce 

J(0) = | dXe- NTr te x '- x -*k(Xy- 1 (9) 

to express the effect on of the fluctuations of the auxiliary matrix (that feels the influence 
of the q — 1 other Potts matrices); writing 

C(Ai) = ^log(J) = <X«> (10) 

we obtain 

2Re - + C(y) = 0,t/G supp a (11) 

This Hartree-Fock formulation of the large N situation is similar to that used in |J to 
find the common distribution of eigenvalues of matrices connected along a Bethe tree; 
here, the tree is finite (with q + 1 sites), and the density is site-dependent. The resolution 
will proceed along similar lines, using results of . 



2 The critical point 

We can extract useful information without solving explicitly (|8|, |Il]) , just by using analytic 
properties of generalized resolvents. Let us introduce 

Kr 'W = j(i^)^ = <(rb)-) (12) 

which we shall also write W x (z), x being the k-th eigenvalue of X, to which is coupled. 
We have § 

1 ^ Wk (z) - Wj(z) l ( d x I 1 x 

n \ — " 'NVT^T^ikk + { J^ )kk (13) 

This formula can be obtained by considering the contracted differential d a bfbc of the 
matrix-valued function: 

/:*-><-!-) 
z — <p 

([□D only expresses the equivariance of / 

f(QXn + ) = Qf{X)Q + 



In the large N limit, ([13]) reduces to 



dxp(x)^ = -w(a)W a (z) + zW a {z) - 1 (14) 

a — x 



Similarly, we define 

1 



Zk{z) ~ j\ z -±dJkk J ^ \\ Z -Xh, 
and we introduce the generalized resolvents 

p(x)W x (z) 



F(z, s) = 1 — dx 



s — X 

aM = LfjyEmM (i 5 ) 

J s-y 

Expliciting F we have 



N —j? S — Xk^^Z — <f) J kk 

l-l/Tr— *— — *— ) (16) 



AM s -X z 

and, as X and change roles from F to G, G(s, 2) = F(z, s). 

For large z, W.(z) is well defined and F(z, •) has a cut along supp p: 

F(z,a±te) = 1 - 4 dx p(y£)Wx(yZ > ±mp(a)WJz) 

J a — x 

= W a (z)[z - w(a) -Re f{a) ±inp(a)} (17) 

in virtue of fll^D- We assume that w and £ can be analytically continued near the supports 
of p, a (independently below and above the cuts: we shall see that the edges are branching 
points) and we set 

u = w + f,v = ( + g (18) 



so that 



F(z,a) = W a (z)(z-u(a)) 

G(z,b) = Z b {z)(z-v(b)) (19) 



and, as these functions go to one at infinity ([Tq) , we can write a dispersion integral for 
their phases 

f dT 1 
F(z, s) = exp - * — log [z - u(t)} 

JSUpp p IVK T — S 

G(z,s) = exp-/ ^-^log[ 2 -,(A)] (20) 

./SUpp a IVK A — S 




Figure 1: In u o r = id, we evaluate u on the sheet on which lies the value taken by r 
after analytic continuation (dashed line). 

with contour integrals going counterclockwise along supp ± ie, this being a valid repre- 
sentation for large s, z. (Proof: for fixed large z, the ratio of both sides is a holomorphic 
function of s, with no cut, going to 1 as s — > oo, thus equal to 1.) We note [a, f3\ = supp p. 

Let us express the symmetry property: G(s,z) = F(z,s): we change the variable in 
the first contour integral from r to u(r) and then integrate by parts to obtain 



To be precise, the former integral involved u± along [a, /3] ±ie; the new expression involves 
the integrals along two arcs, joining u + (a) = w_(a) (for p(a) = 0) to u + {(5) = W-(/3), of 
the corresponding inverse function r: see Figure 1. 

We note that the upper segment is mapped onto the lower arc because Im u + < 0: so, 
the orientation of the image contour depends on the order of u(a), u(f3). 

We expect as few singularities as possible for u and r, because we introduced the 
smallest possible number of coupling constants (should we have taken more, we could 
have produced higher order criticality through a richer singularity structure): we assume 
we can flatten the integration contour to («(«), u{f3)) and analytically continue r without 
encountering any singularity, so that 



We knew that equality for large z, s: under that form, contour integrals are seen to be 
equal (up to 2in, but they vanish at infinity!) for large s and any non-real z. Considering 
the discontinuity on the real z-axis, we obtain for large s: 




(21) 




(22) 



r+{X)-s ] 
r_(A)- a J 



± XAGSUpp a log 



U_(A) - s 



(23) 



where the sign depends on the relative order of u(a), u{(5). We conclude that (u(a), u(f3)) 
covers supp a and that 

r e (A) = f± e (A) , A G supp a (24) 
For instance, if u preserves the order of a, (3: t± = v± oia supp a and 

A = u + [v-(X)] (25) 

for A G supp a and so, for any A by analytic continuation. In any case, attention shall be 
paid to determine the u-sheet on which v (A) has to be evaluated: one has to follow the 
value of t as its argument moves from u(C) towards the real axis. 

This inversion relation u o v = id always holds, in the planar limit, whether we are 
at the critical point or not. However, if v is singular at A then u is at v (A): but at the 
critical point, any function has all its singularities lying at the same place (see section |3]), 
and we expect each density to develop a singularity at one edge of its support; so, these 
singular edges are exchanged under u, v. 

We call 7, 5 the exponents with which p, a vanish at the singular edges: they govern 
the corresponding leading singularities of /, g. For example, if p(a + e) = est e 7 + . . . we 
have 

f(a + C) = regular + M > Q) (-C) 7 + • ■ ■ 

sin 7T7 

where the 7-th power of a positive number is taken positive. Using (§,PJ) we obtain 
similar expansions for u, v, e.g. (singularities at the left edge) 

q- 2 

u±(a + z) = reg + C sin 7r7 [ cos7T7 =1= isimr^z 1 + . . . 

v±(u(a) + z) = reg + C sin n5 [ — cos red =F i sin irS] z s + . . . (26) 

with positive constants Moreover, we expect 1 < 8 < 2 (see 0, @, and section |3|). 

For that singularity to disappear in u o v — id, 7 has to be equal to 5 or to its inverse; 
in the latter case, the linear term in the regular part of v necessarily vanishes. In the 
first case, we write that some linear combination of C, C vanishes and obtain a relation 
between the phases of the bracketed terms of (p6|). In the other case (7$ = 1), the 
requirement that z — {• ■ •)(- 2<S ) 7 a l so gives a relation between these phases. 

Investigating all the possible situations (respective positions of the singular edges, 
relation between 7 and 5) we discover that, for singularities occuring both at the same 
edge: 

either 7 = and cos7r(l — 1/5) = 
or 7 = 5 and cos7r(<5 — 1) = y/q/2 



The case when one singularity occurs at the left edge and the other on the right edge, 
appears to be inconsistent: this is related to our choice of sign of the X.<f> term in (0); 



with that choice X is localized near Y^i'Pii an d we understand that the spectra become 
simultaneously singular at corresponding edges. 

Expressing our results in terms of m — —l/j s tr (parametrization of m, m + 1 unitary 
models) we find 

7 = 5 _1 and = cos7r/(m + 1) 

7 = 5 and = cosn/m (27) 

and recognize the (flat space - regular lattice) critical and tricritical Potts model |J. In 
particular, we recover the unitary (5,6) model as the critical 3-states Potts model. This 
result shows that it is important not to confuse between the exponents of p and sigma, 
otherwise one fails to identify the critical Potts model. 

Before turning to the study of the vicinity of the critical point, a last remark is in 
order: we said that the existence of a singularity for u at A implied a singularity for v 
at u(X). This is true exactly at the critical point, but not in its neighbourhood: out 
of the critical point, densities of eigenvalues vanish as square- roots at the ends of their 
supports (as in the simple one-matrix case 0; see also the scaling functions in the next 
section). Then, u(a + z) is an invertible power series in ^/z, and the condition on the 
inverse function is a zero of its derivative v'(u(a)) = 0. When we approach the critical 
point, that zero disappears at the place where two cuts merge (consider the example of 
y/e 2 — x 2 as e — > 0) and the fractional order of the zero of the density at the edge of 
the support is increased: this will be transparent with the more explicit results of the 
following section. 



3 The scaling regime 

For any values of the temperature and the cosmo logical constant, we can use known results 
about the external field problem (Kazakov and Kostov, in see also |Tl|) to compute 
exactly w in the large N limit: setting a = 1/m, g = g/m 3 , 



V®x + c 1 , 1 /n f pit) dt 1 



... vai + c 1 . f 
wW/a= V_ -j- + 1/2/ 



y/g 2g J yjat + c ^/ax~+~c + \/at + c 

w(x) being the restriction of a holomorphic function with a semi-infinite cut (— oo,— c] 
whose location, at the left of the support of p, is given by 

(for positive g, which we shall assume is the case; for negative g, the cut is located to the 
right of the support of p) 



The eigenvalues of X are in the lowest-energy configuration when 



x 



2 _fdyM + rV^+c_ i / p£ 1 ) ( 30) 



J x — y \ 2(7 7 a/q7 + c v / «a r +~c + + c y 

on the support of p. Changing notations to 



u = \/od + c, tt(u) = p(t) (31) 

and introducing 

the equilibrium equation reads 

2Re G(z) + (2 - q)G(-z) = P(z) (33) 

with 

Of \ ^ ~ ° qa I qa flA\ 

P(z) = -z + — (34) 

Note that 7r is not normalized, but 

while the normalization of p now reads 

u7t(u) du = a/2 (36) 



When q = 2 this equation is easily solved by a dispersion integral and we recover the 
known results about the Ising model on a dynamical lattice, with 

a 2 c = ^^,g/a s = VT0 (37) 

(a c < .5 effectively corresponds to some positive temperature), and the usual scaling law 
for the resolvent (one- loop function). 

For q = or 4, fl33| ) has a Z 2 symmetry (z ^ —z) and its solution is expressed through 
elliptic integrals; we recall that q — > corresponds to the statistics of tree-like polymers 
on a random lattice and has been solved in Jl|. 

In the general case, the critical point is reached when the support of n, [a,b], fuses 
with its mirror image (a — > 0), that is, when the singularity of w reaches the support of 
p. The scaling behaviour of G 

G(acosh6) ~ a^gid) (38) 



is imposed by (|33| ) 

g(9) + g(9 + 2m) + (2 - q)g(9 + in) = 



(39) 



that is 

g(9) = est cosh k9 (40) 

with cos utt = q/2 — 1 which corresponds to the exponent found in the first section for 
the two-critical Potts model. 

So, we found the scaling behaviour of the resolvent of 7r; to obtain the one- loop function 
for a Potts matrix, we shall compute the resolvent of p and use the inversion formula. To 
do so, however, we need to know the position of the support of p which involves c: and 
this parameter is left undetermined when we extract the singular part of G in the scaling 
limit. To assert a definite conclusion about p (and not tt) we need compute c: that we 
did for q = 3, an unsolved model up to now. 

A general way to proceed with (|33l) is to shift G by the polynomial Q that satis- 
fies 2Q(z) + (2 — q)Q(— z) = P(z), to obtain a homogeneous equation for / = G — Q. 
Introducing 9 such that cos# = 1 — q/2 [9 = 27r/3 for q = 3), we define 

Mz) = f(z) + e ±te f(-z) (41) 

These holomorphic functions have a double cut [—b, —a], [a, b], are related by <fi-(z) = 
e~ td (j) + (—z). Through the right cut, 0+ is continued into —</>_; and through the left cut, 
into —e 2ie 4>- according to (|33T). 

When 9 is commensurable to tt, we can raise <f>± to some integer power so as to obtain 
holomorphic functions which realize a double-cover of the sphere (with two cuts). Namely, 
when q = 3, we set ip± = (fi±, so that ip-(z) = ip + (—z) and the analytic continuation of -0+ 
through any cut leads to — ip-. At infinity, ip behaves as z e /a 3 . We obtain a meromorphic 
function on the torus (double cover of the two-slit plane) which can be written 

i)± = a~ 3 (^/(z 2 -a 2 )(z 2 -b 2 )S(z) ± iR(z)) (42) 

with some monic even polynomial 5* of degree 4; and an odd polynomial R with degree 
< 5. 

The vanishing of the sum of the residues of ip' / ip on the torus shows that ip has twelve 
zeroes. These have to be triple zeroes, since = ■0 1//3 has no branching point and we 
conclude that 

S 2 (z)(z 2 - a 2 )(z 2 - b 2 ) + R 2 (z) = (z 2 - \ 2 f{z 2 - p 2 f (43) 

This requirement allows us to solve the Potts model: for given a and g, we look for 
a, b and the coefficients of R, S (3 and 2 parameters); (^3|) gives 4 constraints (for a sixth- 
order polynomial (in z 2 ) to have two triple roots), and 3 other equations are given by the 
asymptotic expansion of / (coefficient of z; eliminate c between the coefficients of z° and 
z~~ l ] coefficient of z~ 2 ). 

2 We thank B.Eynard for his advice 



Precisely at the critical point, we know that a = and p(t c + dt) ~ dt 5 / 6 so that 
tc(u) ~ m 5//3 : inspection of 43 then shows that is a zero of S with multiplicity 4, and 
with multiplicity 5 for R. We thus have identified all coefficients in ^ 



(z 4 ) 2 z 2 (z 2 - b 2 ) + (bz 5 ) 2 = (z 2 ) 3 (z 2 ) 3 (44) 
and we find that 



(45) 



2 _ y/47-3 _ 4 
ttc " 38 ^ 486^/47-3) 

We then solve ( (4~3"D in the vicinity of that point to investigate the scaling limit; to 
obtain a generic perturbation, we can shift g while keeping a = a c . We find the following 
behaviours 

db ~ a 2 , A ~ n ~ a 5/6 with A 2 + /i 2 ~ a 2 (46) 

the roots of R, S scale as a, while the leading coefficient of R varies with a 10 / 3 , just like 
g. And dc = a 2 /2 + . . . 

For the scaling behaviour of / we get 

i r 3 
/(z = acoshy + iO) = j= e~ 2l7r ^ 3 (ia 5 b [sinh y(cosh 4 y — - cosh 2 y + 1/16) 



a 



V3 



— cosh y (cosh 4 y — 5/4 cosh 2 y + 5/16)] )^ 3 
; 2i7r / 3 (ia 5 6 [sinh y(...) + cosh y(. . . )]) 1/3 1 (47) 



where the branch of the third root is fixed by continuation from the behaviour of / at 
infinity, and where the bracketed combinations of hyperbolic lines reduce remarkably to 
— exp(— 5y/16) and exp(5i//16) so that 

f(a cosh 3y + iO) ~ a 5/3 cosh 5{y - 2m/15) (48) 

The explicit limit we found for c allows us to see that p also has a scaling limit expressed 
with hyperbolic lines because u ~ a cosh 3y corresponds effectively to dt ~ a 2 cosh 6y (the 
coefficient 1/2 in (ic = a 2 /2 + ... is really important here) 

Finally, using t>_ ou + = id we find t>_, then the one-loop function for a Potts matrix 
(partition function of a disk with a uniform color along the boundary) 

g(e 5 cosh 5t) = regular + est e 6 cosh 6(£ — in/5) + .. . (49) 

Such scaling laws are universal among unitary models: in || however, a two-matrix 
formulation was used to produce any (m, m + 1) model, with multicritical points analogous 
to Kazakov's one matrix-model (the geometry of the different types of polygons used to 
build the surface -not only triangles- conspires with the Ising spins to produce different 
conformal models) and the identification with the explicit model here discussed is rather 



unclear. We can just argue that in both cases, we compute averages of macroscopic 
insertions of the identity operator. 

To conclude, we would like to consider again the 3-critical point: how could it appear in 
the first section, while we missed it in the above exact computation? To reach tricriticality, 
we need two coupling constants for the statistical model, and not only the temperature: 
it seems that we would obtain a model in the right universality class if we added an TrX 3 
term to the action, to produce a kind of dilute Potts gas; this modification would not 
change the lines of the first section, where the simplest type of criticality was assumed 
with given number of parameters, while it requires a new investigation of the derivation 
given in the last section. That seems interesting to clarify, in order to gain a better 
understanding of the q — > 4 limit, where these two RG fixed points merge into a first- 
order fixed point, this limit being the c — > 1. 
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